Continuous weak measurement of quantum coherent oscillations 
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We consider the problem of continuous quantum measurement of coherent oscillations between two 
quantum states of an individual two-state system. It is shown that the interplay between the 
information acquisition and the backaction dephasing of the oscillations by the detector imposes a 
fundamental limit, equal to 4, on the signal-to-noise ratio of the measurement. The limit is universal, 
e.g., independent of the coupling strength between the detector and system, and results from the 
tendency of quantum measurement to localize the system in one of the measured eigenstates. 



PACS numbers: 73.23.-b,03.65.Bz 

Coherent oscillations between the two states of a quan- 
tum two-state system represent one of the most basic and 
direct manifestations of quantum mechanics and are en- 
countered in practically all areas of physics. The ques- 
tion of how to measure them directly in an individual 
two-state system was formulated Gl-H for the first time in 
the context of quantum dynamics of Josephson junctions, 
where the oscillating variable, the magnetic flux in a su- 
perconducting loop, is macroscopic. A common feature 
of the measurement schemes suggested in this context is 
the use of conventional "projective" measurements that 
localize the flux in one of its eigenstates and suppress the 
oscillations. The time evolution of the oscillations can 
then only be studied if the experiment is repeated many 
times with the same initial conditions, and the informa- 
tion about oscillations is contained in the probability dis- 
tribution of the measurement outcomes. This means that 
the oscillations are effectively studied in an ensemble of 
systems, not in an individual system. Another, practi- 
cal, disadvantage of the projective measurements is the 
need to switch the detector on and off very rapidly, on a 
time scale shorter than the oscillation period, in order to 
allow for free time evolution of the oscillation and sub- 
sequent measurement. Since the oscillation frequency is 
limited from below by several factors, including decoher- 
ence time and temperature, this requirement presents at 
the very least a challenging technical problem. Although 
this problem can sometimes be solved, as demonstrated 
by the recent observation of coherent quantum oscilla- 
tions of charge 0, it is partly responsible for the fact 
that macroscopic quantum oscillations of magnetic flux 
have so far eluded experimental observation. 

The aim of our work is to point out that the problem 
of measurement of quantum coherent oscillations in an 
individual two-state system can be simplified if projec- 
tive measurements are replaced with a weak continuous 
measurement, and to study the quantitative characteris- 
tics of such a measurement scheme. As is emphasized fre- 
quently in the theory of quantum measurements (see e.g., 
|^-^|), the "textbook" projective quantum measurement 
requires the dynamic interaction between the system and 



the detector to be sufficiently strong to establish nearly 
perfect correlation between their states. If the interaction 
is weak, this does not happen, and the measurement pro- 
vides only limited information about the system. Such 
a weak measurement, however, perturbs the system only 
slightly and can be performed continuously. Below we 
consider quantitatively the process of continuous weak 
measurement of quantum coherent oscillations. We cal- 
culate the spectral density of the detector output and 
show that the trade-off between the acquisition of infor- 
mation and dephasing due to the detector backaction on 
the oscillations imposes a fundamental limit, equal to 4, 
on the signal-to-noise ratio of the measurement. In this 
work, we use a more conventional non-selective approach 
to measurement, and all the results can be reproduced 
within the selective description of the measurement pro- 
cess M. 

Although the main conclusions of our work are quite 
general, in what follows we prefer to use the language of a 
particular system: two coupled quantum dots measured 
with a quantum point contact. Quantum point contacts 
were used as electron detectors in [p)|-[i"2] and described 
theoretically in p^-p^j9[|. Coherent electron oscillations 
in coupled dots were observed indirectly in the dc trans- 
port under microwave irradiation ffsl . 

The Hamiltonian of the system (see inset in Fig. la) 

is: 
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The first two terms here describe an electron oscillating 
between the two discrete energy states localized in the 
quantum dots: e is the energy difference between the 
states, —A/2 is their tunnel coupling, and the cr's de- 
note Pauli matrices. The operators represent point- 
contact electrons in the two scattering states i = 1, 2 (in- 
cident from the two contact electrodes) with momentum 
k. The coupling a z U/2 is due to an additional scatter- 
ing potential ±U(x)/2 created in the point contact by 
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the electron occupying one or the other dot. The point 
contact is biased with a dc voltage V, so that changes 
in the scattering potential lead to changes in the current 
/ through the contact. We take eV to be much smaller 
than both the Fermi energy in the point contact and the 
inverse traversal time of the contact. This allows us to 
linearize the energy spectrum of the point-contact elec- 
trons: £fc = Vpk, where vp is the Fermi velocity, and 
neglect the momentum dependence of the coupling ma- 
trix elements Uij = J dx ip* (x)U (x)ipj (x) of the pertur- 
bation U in the basis of the two scattering states ipi(x). 
We also assume that the Uij are sufficiently small for 
the point contact to operate as a linear detector, and 
treat the contact's response to electron in the dots in the 
linear-response approximation. 

Quantum oscillations of electron between the dots cre- 
ate an oscillating component of the current / through the 
point contact. Since the phase of the oscillation diffuses 
under the backaction of the shot noise of the point con- 
tact, the oscillations are best characterized by their spec- 
tral density. To find the spectral density of the current I 
we choose the origin of the coordinate x along the contact 
in such a way that the unperturbed scattering potential 
is effectively symmetric, i.e., the reflection amplitudes for 
both scattering states are the same. Then, the current 
operator calculated at a point x in the asymptotic region 
of the scattering states is: 

I= evF Y^[D(a{ k a lp ~ a\ k a 2p ) + 
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where D and R = 1 — D are the transmission and reflec- 
tion probabilities of the point contact, L is a normaliza- 
tion length, and the variation of the momentum near the 
Fermi points (i.e., the difference between k and p) was 
neglected everywhere besides the phase factor in the sec- 
ond term. The reason for keeping this factor will become 
clear later. 

In the linear-response regime, the current response of 
the point contact is driven by the part of the perturbation 
U causing transitions between the two scattering states 
"01,2 ■ Considering the effect of this perturbation on the 
stationary (symmetric and antisymmetric) combinations 
of the scattering states, one can show that the real part of 
the transition matrix element U\i is related to the change 
5D of the transmission probability of the contact: 
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The imaginary part of U12, expressed through a dimen- 
sionless parameter u in eq. (||), does not affect the current 
I. Qualitatively, it characterizes the degree of asymmetry 
in the coupling of the quantum dots to the point contact; 
u = if the perturbation potential U{x) is applied sym- 



metrically with respect to the main scattering potential 
of the point contact. 

When the point contact is used as a detector in a quan- 
tum measurement, the current / plays the role of the 
measurement output and should behave classically. This 
condition requires the spectral density of I to be much 
larger than the spectral density of the zero-point fluctua- 
tions in the relevant frequency range. It is satisfied when 
the voltage V across the point contact, which determines 
the magnitude and the threshold frequency of the shot 
noise of J, is sufficiently large, eV ^> e, A. For the point 
contact to be an effective detector, eV should also be 
much larger than the temperature T . In this regime, it 
is straightforward to find the correlation functions of the 
perturbation U and the current / in the zeroth order in 
V from eqs. (|), (§), and (§): 



e 2 V , 



(4) 



(U(t)I(t + r)) = (iSD + u) 5{t - rf) . (5) 

The spectral density of I at frequencies much smaller 
than eV is dominated by the regular shot noise, and the 
current correlation function is Kj ^ (t) = (I(t+T)I(t))o = 
e(I)R5(r), where (I) = e 2 VD/n. The time delay 77 = 
I^I/vf in eq. (0) comes from the phase factor e~ l ( k -p)\ x \ 
kept in eq. (0), and is infinitesimally small for small 
traversal time of the contact. It is nevertheless impor- 
tant for resolving the ambiguity in averages involving the 
time ordering of I and U that are needed for the calcu- 
lation of the current response: i J dt' (T {I (t)U (i')}) = 
e 2 V{5D + iu)/ir. 

Expression for the current correlation function Kj(t) 
in the interaction representation with respect to U is: 

K I (T)=Tl-{p{t)I(t)S\t + T,t)I(t + T)S{t + T,t)}, (6) 

where p(t) is the total density matrix of the point contact 
and quantum dots at time t, the trace is taken over both 
systems, and S(t + r,t) = Texp{(— i/2h) f* +T dt'a z U} 
is the time evolution operator. Taking the trace over the 
electron states in the point contact in eq. (^|) with the 
help of the correlation functions (Q), we get 



K I (r) = Kf\r) 



(SI) 2 



(7) 



The average ( • • • ) in eq. (fjj) is taken over the two states 
of the quantum dots with the stationary dot density ma- 
trix p established as a result of the interaction with the 
point contact and averaged over its dynamics. The cur- 
rent change 81 = e 2 (5D)V/ir is the current response to 
electron oscillations between the dots, and cr z (r) now de- 
notes the full time evolution of a z , driven both by the dot 
Hamiltonian and the interaction U with the point con- 
tact. Qualitatively, eq. (Q) shows that the current corre- 
lation function directly reflects the correlation function 
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of the electron position in the dots given by the operator 
er z . 

The time dependence of the operator <t 2 (t) in eq. (0) 
is obtained by tracing out the point contact degrees of 
freedom in eq. (^) with the help of the U-U correlation 
function (^) . In this way we get the standard set of equa- 
tions for the matrix elements ov, of a z (r) : 



&u = AImcr 12 , &12 = (is — r)cri2 — iAoii 
and 022 = —(fit- The rate 



T = eV 
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describes backaction dephasing of the coherent electron 
oscillations between the dots by the point contact. Equa- 
tion (^|) shows that the dephasing rate reaches a min- 
imum in the case of symmetric dot-contact coupling 
(u = 0). In this case, the rate of dephasing by a point 
contact has been found in fl3|-|l5[|. Increased dephasing 
in the case of asymmetric dot-contact coupling was dis- 
cussed qualitatively in [fl7| and studied experimentally 
in . Since the decrease of T with decreasing asymme- 
try u does not affect the current response of the point 
contact, symmetric coupling corresponds to an optimum 
in its operation as a detector. In this regime, the point 
contact represents an ideal quantum detector in a sense 
that the minimum value of the dephasing rate (|^) is de- 
termined purely by the rate of information acquisition 
about the state of the quantum dots and can be written 
as r = (5I) 2 /4So, where So = 2e(I)R is the spectral den- 
sity of the current shot noise of the point contact || . This 
part of the dephasing is fundamentally unavoidable and 
reflects the tendency of quantum measurement to local- 
ize the measured system in one of the eigenstates of the 
measured observable, in our case, the electron position 

The dot density matrix p satisfies the same set of equa- 
tions (|8|), except for the normalization, P11+P22 = 1, and 
its stationary value is p = 1/2. Solving eqs. (|^) with the 
initial condition a z (0) = a z and averaging <j z g z (t) over 
p = 1/2 we find the correlation function (0) and the spec- 
tral density Sjfw) = 2 Ki(r)e lu}r dr for e = 0: 



Si{lo) = So 



Til 2 (SI) 2 



(tu 2 -n 2 ) 2 + v 2 w 2 



(10) 



In the case of biased dots with e ^ 0, it is convenient 
to calculate the spectrum numerically from eq. (||). The 
spectrum in this case is plotted in Fig. 1 for several values 
of e and the dephasing rate T. For weak dephasing, T <C 
A, the spectrum consists of a zero-frequency Lorentzian 
that vanishes at e = and grows with increasing |e|, and 
a peak at the oscillation frequency ft = (A 2 + e 2 ) 1 / 2 . 
Although the width of the oscillation peak is T and can 
be small for sufficiently weak dot-contact coupling, its 



height cannot be arbitrarily large in comparison to the 
background noise spectral density Sq. At e = 0, when 
the amplitude of the oscillations is maximum, the peak 
height is S max = (SI) 2 /T. Even in this case, the ratio of 
the peak height to the background is limited: 



A(8Df 



So (SDy 



< 4. 



(11) 



This limitation is universal, e.g., independent of the cou- 
pling strength between the dots and the point contact, 
and reflects quantitatively the interplay between mea- 
surement of the quantum coherent oscillations and their 
backaction dephasing. The fact that the height of the 
spectral line of the oscillations can not be much larger 
than the noise background means that, in the time do- 
main, the oscillations are drowned in the shot noise. 
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FIG. 1. The diagram (inset in (a)) of the coherent electron 
oscillations between the two discrete energy states in coupled 
quantum dots measured with a point contact. The oscillations 
are detected through modulation of the current I(t) in the 
point contact biased with a voltage V. Plotted curves show 
the spectral density Si(u>) of the current I in the case of 
symmetric coupling between the point contact and the dots 
for several values of (a) the energy bias e between the dots 
reflected in the oscillation frequency Q = (A 2 + e 2 ) 1//2 , and 
(b) the rate F ([]) of the measurement-induced dephasing. 
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The total intensity of the oscillation line in the spec- 
trum: 

. dw (SI) 2 



[Si(w) -S ] 
o 27r 4 



(12) 



does depend on the strength of coupling to the point con- 
tact, increasing as the coupling becomes stronger. An 
interesting feature of eq. ( |l2|) is that it stresses the im- 
possibility of a simple classical interpretation of the quan- 
tum coherent oscillations, since the intensity of harmonic 
classical oscillations of the same amplitude SI/2 would 
be two times smaller, and no classical signal of this am- 
plitude could produce the oscillation line with intensity 

When the backaction dephasing rate T increases, the 
oscillation line broadens towards the lower frequencies, 
and eventually turns into the growing spectral peak at 
zero frequency that reflects the incoherent electron jumps 
between the two dots. At large T, when the coherent os- 
cillations arc suppressed, the rate of incoherent tunneling 
decreases with increasing T. For instance, at T ^> f2, 
the tunneling rate is 7 = A 2 /2T , and the spectral 
density of the point contact response has the standard 
Lorentzian form Si(u) - S = 2-f(SI) 2 / (4-f 2 + uj 2 ). Sup- 
pression of the tunneling rate 7 with increasing dephas- 
ing rate T is an example of the generic "Quantum Zeno 
Effect" in which quantum measurement suppresses the 
decay rate of a metastable state. In the context of search 
for the macroscopic quantum coherent oscillations, the 
Lorentzian spectral density has been observed and used 
for measuring the tunneling rate of incoherent quantum 
flux tunneling in SQUIDs |l9| . 

The maximum signal-to-noise ratio S max /So ([ll]) is 
attained if the fundamental backaction of the detector 
is the only dephasing mechanism of the coherent oscil- 
lations. In the case of measurement with a point con- 
tact, the fundamental measurement-induced dephasing 
considered above is created by the backscattering part 
U12 (0) of the dot-contact interaction that dominates at 
large bias voltages V. The forward scattering Un, U22 
does not affect the current / in the contact but creates a 
weak additional dephasing and energy-relaxation mech- 
anism for the oscillations. We now want to discuss the 
effect of such a weak relaxation on the spectral density 
of the oscillations noticeable if the backaction dephasing 
is also weak, T A. 

The inclusion of the additional weak relaxation does 
not modify the calculations that lead to eq. (]?]), 
apart from a trivial modification that now the aver- 
age a z is non-vanishing, and the current correlation 
function should be calculated as Kj(t) = Kj°\t) + 
(SI/2) 2 [(l/2)(a z a z (T) + a z (r)a z ) ~ (cr 2 ) 2 ]. For weak cou- 
pling, it is convenient to find the time evolution of c 2 (t) 
in the basis of eigenstates of the two-state Hamiltonian 
— (ea z + Aa x )/2. Solving the Heisenberg equation of mo- 
tion up to the second order in the dot-contact coupling, 



and tracing out the contact degrees of freedom, we get a 
set of equations for the evolution of the matrix elements 
Sij of o z (t) in the eigenstates basis: 

s 3] {t) = T e [^ - coth{^} Sji ] + (-l) J ^(sn - S22) , 



L n L 2T J ni 

si2(t) = {is - To)si2 



2fl 2 



(13) 



with the initial conditions sn = — S22 = c/^, and S12 = 
— A/f2. The characteristic energy-relaxation rate in eq. 
@ is T e = vA 2 /fl, where v = (l/7r)([/ 1 2 1 + t/| 2 )(L/u F ) 2 , 
and the total dephasing rate is 

T = [w(A 2 ficoth(fi/2T) + 4e 2 T) + r(2e 2 + A 2 )]/2ft 2 . 

The dot density matrix r in the eigenstates basis satisfies 
similar equations, and the stationary values of its matrix 
elements are rn — (T e + r t )/2r t and ryi — 0, where 
T t = r e coth(fi/2T) + TA 2 /tt 2 . From these relations 
and eqs. ( |l3|) we find the spectral density: 



Jfi_f££)2i_£* 
1 [ T t n uj 2 + 



A 2 



E 



(w ± 



r 2 
1 



(14) 



As before, the spectral density consists of a zero- 
frequency Lorentzian and peaks at ±fl of width Tq that 
represent the coherent electron oscillations. Energy re- 
laxation with characteristic rate T e broadens the oscilla- 
tion peak and reduces its height S max , so that the relative 
magnitude of the peak, S max / Sq decreases in comparison 
with its value without relaxation. 

In summary, we have considered a continuous weak 
quantum measurement by a point contact of quantum co- 
herent oscillations in a two-state system, and calculated 
the spectral density of the output signal of the measure- 
ment. It has been shown that the backaction dephasing 
introduced into the oscillation dynamics by the measure- 
ment imposes the fundamental limit on its signal-to-noise 
ratio. We also calculated the effect of energy relaxation 
on the output spectrum. 
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